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ABSTRACT 



Mass loss appears to be a common phenomenon among disk-accreting astrophys- 
ical systems. An outflow emanating from an accretion disk can act as a sink for mass, 
angular momentum and energy and can therefore alter the dissipation rates and ef- 
fective temperatures across the disk. Here, the radial distributions of dissipation rate 
and effective temperature across a Keplerian, steady-state, mass-losing accretion disk 
are derived, using a simple, parametric approach that is sufficiently general to be 
qs ' applicable to many types of dynamical disk wind models. 

Effective temperature distributions for mass-losing accretion disks in cataclysmic 
variables are shown explicitly, with parameters chosen to describe both radiation- 
driven and centrifugally-driven outflows. For realistic wind mass-loss rates of a few 
percent, only centrifugally-driven outflows - particularly those in which mass loss is 
concentrated in the inner disk - are likely to alter the disk's effective temperature dis- 
tribution significantly. Accretion disks that drive such outflows could produce spectra 
and eclipse light curves that are noticeably different from those produced by standard, 
O ' conservative disks. 
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1 INTRODUCTION 

Accretion disks are an important ingredient in our current understanding of many astrophysical systems on all scales. Examples 
of presumed disk accretors include young stars, compact objects in close binary systems and active galactic nuclei (AGN) and 
quasars (QSOs). There is also evidence that the process of mass accretion via a disk is often - and perhaps always - associated 
with mass-loss from the disk in the form of a wind or a jet. For instance, the broad- line region in QSOs, the radio jets in 
AGN, the jets and bipolar outflows from young stars and the fast winds observed in disk-accreting cataclysmic variables all 
provide compelling evidence for such a disk- wind connection [see Livio (1999) for a recent review]. 

Since the outflow provides an additional sink for mass, angular momentum and energy, the spectrum emitted by a mass- 
losing accretion disk will in general be different from that produced by a "conservative" one. Intuitively, mass loss may be 
expected to reduce the effective temperatures across the entire disk, but especially in those regions from which matter is 
preferentially expelled. Based on this expectation, it has been suggested, for example, that mass loss from the inner disk 
may be responsible for the apparently flatter-than-expected radial brightness temperature distributions inferred for some 
disk-accreting nova-like CVs (Rutten, vanParadijs & Tinbergen 1992). 

Despite this, few attempts have been made to generalize steady-state accretion disk theory to include the family of mass- 
losing disks. As a result, the standard T oc R~ ' law that is valid for conservative disks is almost invariably used in practice 
(i.e. in modeling and in comparisons between theory and observations), even in the context of disk-accreting systems in which 
mass loss could alter the disk's effective temperature distribution significantly. Where modifications to the standard picture 
have been attempted, these have usually been ad hoc (e.g. Shlosman, Vitello & Mauche 1996). An exception is the work by 
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2 C. Knigge 

Piran (1977), which provides analytic expressions (in integral form) for the energy generation rate in a mass- losing accretion 
disk. However, the bulk of Piran's paper deals specifically with thermally-driven, evaporative disk winds. The only other type 
of mass-loss considered explicitly by him is a "delta-function outflow" from a particular radius on the disk. 

The purpose of the present work is twofold. First, simple, analytical expressions are derived for the radial distributions 
of dissipation rate (Section |2|) and effective temperature (Section |3|) across a steady-state accretion disk that is losing mass 
and angular momentum to an outflow from its surface. The derivation is designed to closely parallel standard, steady-state 
accretion disk theory, as described, for example, by Pringle (1981). The parametric model adopted to describe the mass loss is 
simple, yet sufficiently general to be applicable to many types of dynamical disk wind models, including both radiation-driven 
outflows (e.g. Murray & Chiang 1996; Proga, Stone & Drew 1998; Feldmaier & Shlosman 1999; Feldmaier, Shlosman & Vitello 
1999) and centrifugally-driven, magneto-hydrodynamic (MHD) disk winds (e.g. Blandford & Payne 1982; Cannizzo & Pudritz 
1988; Emmering, Blandford & Shlosman 1992; Pelletier & Pudritz 1992). 

Second, the formalism is applied to accretion disks in cataclysmic variable stars (CVs; Section rcS(). The mechanism respon- 
sible for the mass loss from these systems is currently still unknown, so both radiation-driven (Section p.l|) and centrifugally- 



driven outflows (Section 3.2) are considered as possible culprits. It is found that if CV winds are of the latter type, they could 



leave an observable imprint on the disk's temperature profile. 



2 THE VISCOUS DISSIPATION RATE IN A STEADY-STATE, MASS-LOSING DISK 

It is convenient to start by defining the cumulative mass-loss rate from the disk as 

M W (R)= rh w {R') A-kR' dR . (1) 

Here i?* denotes the radius at the inner edge of the disk (R* ~ Rwd in non-magnetic CVs, where Rwd corresponds to the 
radius of the white dwarf [WD] at the centre of the disk) and m m (R) is the mass-loss rate per unit area from each disk face. 
With this definition the total mass-loss rate is given by 

M Wi tot a l = M m (Rdisk), (2) 

where Rdisk is the radius of the accretion disk. The radial rate of change of M W (R) is just 

dM w (R) 

— ' = 4nRm m (R). (3) 

The equations of mass and angular momentum conservation for a thin disk subject to mass and angular momentum loss 
can now be derived in close analogy to the usual case by considering a particular annulus on the disk with inner radius R and 
outer radius R + AR (Pringle 1981; Frank, King & Raine 1992). 

The mass of this annulus is 2irRARE - where £ is the surface density of the disk - and the time rate of change of this 
is equal to the rate at which mass flows into the annulus from its radial borders minus the rate at which it loses mass to the 
wind. In the limit AR — > 0, mass conservation then gives 

dt R 8R ( R> 2nR dR ' ( ' 

where Vr is the radial inflow velocity of material in the disk (Vr < 0). 

In order to write down the angular momentum equation, it is necessary to specify how, and how much, angular momentum 
is lost to the wind. Here, it will be assumed that matter ejected at radius R on the disk carries away specific angular momentum 
(IR) Q., where fi is the angular velocity in the disk at R. Thus / = corresponds to a non-rotating disk wind and I — 1 to 
outflowing material which retains the specific angular momentum it had at the point of ejection. Models in which the lever-arm 
/ > 1 effectively describe an outflow with bead-on-a-wire rotation out to a radius Rout = IR, as envisaged for centrifugally- 
driven magnetic disk winds. Models such as this can remove a lot of angular momentum from the disk. 

With the angular momentum carried by the outflowing material thus fixed, conservation of angular momentum gives 

—dT~ ~ -RM {R S%n) + R8R {R vS dR> ~ 2nR dR " (5) 

The first two terms on the right hand side of this expression describe the inflow of angular momentum through the boundaries 
of the annulus and the effects of viscous torques due to shear (v here is the effective kinematic viscosity). The third term 
allows for the angular momentum sink provided by the outflow. Expressions closely analogous to Equations M) and (Eh were 
first written down by Bath, Edwards & Mantle (1983), who examined the process of stream penetration, i.e. the effects on 
the disk of a radially varying mass input with constant specific angular momentum. 
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Steady- state, mass-losing accretion disks 3 

To keep the discussion simple, the disk will now be assumed to be Keplerian and in a steady-state. The first of these 
assumptions gives Q(R) = (GMt/R 3 ) 1 ' 2 , where M* is the mass of the accretor, whereas the second implies that J^ = in (^ 
and (|J). The conservation equations can then be integrated, which for (||) gives 



2nRT,V R - M W (R) - constant. 



(6) 



The first term on the left hand side of this expression is clearly the accretion rate at radius R on the disk, M acc (7?). The 
constant of integration in (a) can be found by considering the boundary condition at R — R, , which shows that it is simply 
the rate of accretion onto the central star, i.e. 



MaUR) = M acc (R*) + M W (R). 



(7) 



The angular momentum integral is somewhat more difficult to compute, since as a result of the factor of R Q in the new 
outflow sink term, the right hand side of Equation (pi) is no longer a perfect differential. That term, at least, must therefore be 
integrated explicitly, which requires M m (R) or equivalently rh m (R) to be specified. To keep things tractable, a simple power 
law form will be used here, 



rh w (R) = KR i , 



(8) 



with the constant K being fixed by normalizing to M Wtt otai (which is used as a free parameter). Note that, in principle, the 
local mass-loss rate may itself depend on the local effective temperature; in this case an iterative procedure would be needed 
to deal with the resulting non-linearity. 

Equation (|5f) can now be integrated straightforwardly. Using the usual boundary condition, — i*' = (which really 
applies at or near the boundary layer, rather than at the stellar surface itself; e.g. Frank, King & Raine 1985), and the 
identification — 2-kR'EVr = M aC c(R) made above for the mass conservation integral, this yields 



Macc(R) ~ Macc(R*) (^f ) + Mw,total ( ^f ) 
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(9) 



Now the rate at which energy is generated by viscous shear in the disk (per unit area on each disk face) is 



dRj V 87ri? 3 J 

(e.g. Frank, King & Raine 1985). Together with Equation (H), this yields 



D(R) = { 



(W-) M acc {R) - M acc (R») (^f) 1/2 + M Wttotal (^) 1/2 
(ffr-) lM acc (R) - M acc (R«) (£ ) 1/2 + M w , total (^) 1/2 [2/ 2 
(f^)(M a c C ( J R)-Mac C (^)(^) 1/2 + M» )to ^(^) 1/2 
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,1/2 



C^-2,-5/2 



(11) 



;*) 



£ = "5/2. 



for a Keplerian disk. This expression correctly reduces to the usual result for the case of no mass loss upon setting M Wt total — 
and noting that M a cc{R) = Macc(Rdisk) — const in a conservative disk. The explicit form of Equation |ll| for a conservative 
disk is 



Dfio wind\Fi') 



3GM,M ac 
8nR 3 



R*\ 1/2 



(f) 



(12) 



It is convenient to rewrite Equation [ill in such a way that the difference between it and the conservative case, Equation Q2J, 
becomes more explicit. After some algebraic manipulation, one finds that D(R) can be expressed as 



D(R) = D no wmd {R) - 3GM,M totai x 

oivR^ 
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where 



X(R) = < 



l-(t) 
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For reference, it is noted without explicit derivation that if m(R) is a delta-function centered on Rw, then D(R) can still 
be written in the form of Equation [13l but with X(R) now given by (c.f. Equation 2.6 in Piran 1977; note that a power 1/2 
is missing from the second term in the 7? > Rw part of his equation) 



'Rw\ 1 / 2 



1/2 



X(R) = 



i*m h ' -m' *>r 



i-(*) 



1/2 



(15) 



R < R W - 



3 THE RADIAL EFFECTIVE TEMPERATURE PROFILE OF A STEADY-STATE, MASS-LOSING 
ACCRETION DISK 

Given the rate of viscous dissipation in the steady-state disk, the radial effective temperature profile T(R) is usually found 
by assuming that the rate at which accretion energy is released by dissipation is balanced by radiative losses, i.e. 



D{R) r 



°Ttff{R) 



(16) 



if the disk is optically thick. However, if a wind is driven from the surface of the disk, some of the dissipated accretion energy 
may not be radiated away to infinity, but may go instead into powering the wind. In particular, two energy sinks need to be 
considered: (i) the binding energy that must be overcome to allow matter to escape at all; (ii) the kinetic energy carried by 
the flow at infinity, [j 

Strictly speaking, a self-consistent dynamical model of the full disk + wind system is needed to properly describe the 
process of converting the energy dissipated in the disk into wind energy. However, in the spirit of the current treatment, whose 
point is to elucidate the effects of mass loss on the disk's energy balance more generally, it makes sense to simply treat wind 
energy losses, l w (R), as a local cooling term, i.e. 



D(R) = aT* ff (R)+l w (R). 



(17) 



Local approximations such as this are crude, since they do not allow for the possibility that the dissipated accretion 
energy may be redistributed radially before powering the wind. It is also worth noting that, even though Equation |l7] can be 
made energetically correct by construction (see below), it is quite possible for l w (R) to have a strong wavelength dependence. 
For example, in a line-driven wind, energy losses from the disk's radiation field to the wind may be concentrated in a few 
strong spectral lines. In this case, estimates of the effective disk temperature derived from the observed continuum would 
yield values corresponding to l w (R) ~ in Equation |l7|. 

Ignoring these subtleties for the moment, some general statements can be made about the conditions under which a wind 
loss term is required in the energy balance equation at all. First, the rotational energy that is extracted from the disk by 
the outflow is insufficient to overcome the binding energy of the outflowing matter unless I 2 > 3/2. This can be shown, for 
example, by integrating D(R) over R for the case where rn(R) is a delta function (Equations h3J and h5|) and comparing the 
result to the case of a conservative disk. This yields a "luminosity depletion" of (c.f. Equations 2.7 and 2.8 in Piran 1977) 



AL- 



GMM U 



2R, 



(3/2 -f) 



GMM, 



w, total 



Rw 



(18) 



* By not explicitly accounting for the thermal energy content of the wind, we are effectively considering only models in which all of the 
wind's thermal energy is eventually converted into bulk kinetic energy (i.e. models in which T mind ->0asr-i oo). Close to the disk, 
the thermal wind energy must in any case be small compared to the binding energy unless the mass loss is in fact thermally driven. 
This requires temperatures on the order of T esc ~ 3 X 10 8 (M* / Mq)(R/ 10 9 cm)~ ,J K near the base of the flow. Such outflows have been 
considered by Piran (1977) and by Czerny & King (1989ab). 
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Steady- state, mass-losing accretion disks 5 

If the wind material had never been accreted at all, the depletion would have been equal to the first term in Equation [18|. 
Thus, unless I 2 > 3/2, the disk still gains energy from the matter that is ultimately lost to the outflow; for I = 1, this energy 
is just the binding energy of the wind material. This result should not come as a surprise: for 1 = 1, the outflow clearly exerts 
no net torque on the accretion disk, and thus no energy is extracted from the disk to actually power the wind. The accretion 
disk therefore dissipates all of gravitational binding energy it extracted from the wind material as it accreted from infinity to 
Rw prior to its ejection. For I < 1, even more energy can be dissipated by the disk, since the outflowing material carries with 
it less rotational kinetic energy than it had at the point of ejection. Specifically, for I — 0, the disk gains both the rotational 
kinetic energy of the wind material (GMM w /2Rw) and its potential energy (GMM W /Rw)- 

Thus unless an external energy source is invoked to help power the wind, l w must, in general, contain at least an effective 
binding energy term for models with I 2 < 3/2. After all, the outflowing material cannot reach infinity unless at least its 
binding energy (including any rotational energy it "injected" into the disk) is returned to it. Equation hq also shows, however, 
that models with I 2 > 3/2 do extract sufficient energy from the disk to overcome the binding energy. In fact, if I 2 » 3/2, 
the energy that is extracted from the disk is sufficient to power the wind to potentially very high terminal velocities of 
v oo = -v/2(P — 3/2)GM*/Rw — \/2lGM,/Rw- This illustrates how centrifugally-driven accretion disk winds are ultimately 
powered by accretion energy and shows that explicit binding and kinetic energy terms need not be added to the right-hand 
side of Equation hq when modeling disks that drive outflows of this type. 

With these general points in mind, the wind loss term, l w , can be conveniently parameterized as follows. The rate at 
which energy must be supplied (per unit area) to wind material in order for it to overcome its binding energy and escape 
from the disk is 

h{R) = ^m(R)v 2 K (R), (19) 

where vk(R) denotes the Keplerian velocity. Here, no distinction has been made yet between models with different lever arms, 
which can have very different effective binding energies; this will be accounted for below by means of an efficiency factor. 
Similarly, kinetic energy has to be supplied to the wind at a rate of 

l k (R) = ^m(R>l(R) = \m(R)f 2 v 2 K {R) (20) 

per unit area, where the natural scaling v x (R) oc v f , sc (R) oc vk(R) has been assumed, with / a constant. Here, v esc (R) 
denotes the local escape velocity at a given streamline foot point. 

Depending on the length of the lever arm I and on whether there is an external energy source that helps to power the 
wind, both It and Ik may contribute fully or in part to l w . Allowing for this, l m can be written as 

lw(R) = Vbh(R) + vMR) = \(nb + Vkf)m w (R)v 2 K (R). (21) 

where rjt and r)k are efficiency factors. If there are no external energy sources, 

3 - 2l 2 I 2 < 3/2 
,„. = { (22) 

I 2 > 3/2 



and 



I 2 < 3/2 



% = < „ „ . (23) 



(f -3/a)vi., _ 1 2(i 2 -3/2) 



V > 3/2. 



Equation NT] with l w (R) given by Equation Ell is the sought-after expression for the radial effective temperature distribution 
across a mass- losing accretion disk. On substituting dill) or dT3) into (pj), it is easy to solve for T e ff(R) corresponding to any 
consistent set of parameters. More explicitly, the result is 



4 3GM.M ne 

Feff{R > - 8nR* 



, ,-R.\ 1/2 



(f) 



3GM*M W total ■ v ,r,\ 1 / r2x • / n \ 2 , ni , r, A \ 

g^3 X{R) --(» + Vkf )m w (R)v K (R), (24) 



where X(R) is given by Equation (114). The physical meaning of the terms on the right-hand side of this expression is 
straightforward: the first term is the dissipation rate in the absence of any outflow; the second term accounts for the reduction 
in dissipation that results from the loss of mass and angular momentum from the disk to the wind; the third term represents 
the cooling that results if the energy dissipated by the disk must also supply some or all of the wind's binding and kinetic 
energies. Within the second term, in the expression for X(R), terms that involve (do not involve) the lever arm I arise because 
the wind acts as a sink for angular momentum (mass). 
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6 C. Knigge 

4 SUMMARY OF THE MODEL PARAMETERS 

A total of six parameters are needed to specify how the presence of the accretion disk wind affects the disk's radial effective 
temperature distribution: 

(f ) The length of the rotational lever arm I. This variable permits the most meaningful division of the available parameter 
space and, in principle, allows three types of accretion disk winds to be identified: 

(i) I < 1: This corresponds to the family of disk winds that carry away less angular momentum than possessed by the 
wind material before it left the disk surface. Thus each disk annulus loses a larger fraction of its mass than of its 
angular momentum to the outflow. Even though such disk winds provide a net angular momentum sink (unless 
I — 0), they actually increase the specific angular momentum of the remaining disk material. This family of outflows 
is unlikely to be physical and will not be considered further below. 

(ii) I — 1: This corresponds to the family of specific angular momentum-conserving disk winds. This prescription should 
be most appropriate for radiation-driven outflows (e.g. Murray & Chiang f 996; Proga, Stone & Drew f 998; Feldmaier 
& Shlosman f999; Feldmaier, Shlosman & Vitello f999). 

(iii) I > 1: This corresponds to the family of accretion disk winds that remove both specific and net angular momentum 
from the disk. This prescription is appropriate for centrifugally-driven MHD winds (e.g. Blandford & Payne f982; 
Cannizzo & Pudritz 1988; Emmering, Blandford & Shlosman 1992; Pelletier & Pudritz 1992). In these, the length 
of the lever arm corresponds to the Alfen radius, Ra, out to which bead-on-a-wire type rotation on magnetic field is 
(effectively) enforced.FI Thus I = Ra/R in this case. Note that, in general, thermally-driven outflows also belong to 
this class (Piran 1977). 

(2) The total wind mass-loss rate, M Wi totai- 

(3) The radial mass loss power law index £. 

(4) The velocity parameter / = v^/vk- 

(5) The fraction rji, of the wind's binding energy that is provided by dissipated accretion energy. Note that for I < 1, rjt, can 
be greater than unity; see Section H. 

(6) The fraction r\k of the wind's kinetic energy at infinity that is provided by dissipated accretion energy. 



5 CONSISTENCY OF PARAMETER COMBINATIONS 

The purely parametric approach taken above to obtain the dissipation rate and effective temperature distributions across a 
steady-state, mass-losing accretion disk does not ensure that a self-consistent solution exists for any given set of specified 
parameters. This is easily seen, for example, by noting that nothing prevents us from specifying the mass-loss rate to be 
larger than the accretion rate, even though it is clear that a steady-state solution cannot exist in this case. Similarly, real disk 
winds may not give rise to power law mass loss distributions. Once again, a full dynamical model of the disk+wind system 
would be needed to avoid such pitfalls. The parameterized modeling approach adopted here is nevertheless useful because it 
provides a simple and general means for estimating how, and how much, different types of "real" outflows would affect a disk's 
effective temperature distribution. Such estimates may be important as input to other types of kinematic modeling efforts 
(e.g. Shlosman, Vitello & Mauche 1996; Knigge & Drew 1997) and can aid in the interpretation of T e ff(R) distributions 
derived from observations (e.g. Home 1985; Rutten, vanParadijs & Tinbergen 1992). 

One consistency check that can be performed even within the present parametric approach is to identify energetically 
inconsistent parameter sets. These reveal themselves by giving rise to formally negative dissipation rates and effective tem- 
peratures in some or all disk regions, and, strictly speaking, all such solutions should be rejected outright. However, some 
parameter sets may only just fail to achieve energetic consistency, for instance as a result of the turnover in the (standard) 
dissipation rate distribution in the innermost disk regions, which the adopted power law mass loss distribution fails to follow, 
ft is therefore useful to define a quantitative measure of the degree of energetic inconsistency exhibited by a given solution 
and to consider a somewhat less restrictive rejection criterion based on this. 

One such measure is the ratio, 1Z, of the "negative luminosity" contained in any disk regions in which the dissipation rate 
and/or effective temperature is formally negative to the positive luminosity that remains in all other disk regions. The absolute 
value of this ratio is the fraction of the remaining luminosity that would have to be diverted to the negative temperature 
regions in order to make the solution consistent there. We may then choose to accept solutions for which \1Z\ 7^ 0, provided 
that \R\ « 1. The plausibility argument behind this more lenient rejection criterion is that only a minor adjustment of the 

T In reality, wind material does not co-rotate all the way out to Ra- However, the total specific angular momentum carried away by 
such an outflow (including magnetic torques) is exactly Q,R 2 A . So, as far as the angular momentum - and the disk - are concerned, the 
flow effectively co-rotates out to Ra (e.g. Spruit 1996). 
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Steady- state, mass-losing accretion disks 7 

disk structure would probably be needed for such a solution to become consistent and that the (non-negative) dissipation rates 
and effective temperatures across the disk would only have to change slightly to support this adjustment. For all parameter 
combinations examined in the following section \TZ\ < 0.03, and for most combinations \TZ\ = exactly. 



6 MASS-LOSING ACCRETION DISKS IN NON-MAGNETIC CATACLYSMIC VARIABLES 

In this section, the formalism developed above will be used to consider the effect of mass loss in the form of radiation- 
driven (Section |6.l| ) and centrifugally-driven winds (Section p.2J ) on the effective temperature distribution of accretion disks 
in cataclysmic variable stars (CVs). These two types of outflows show the most promise for explaining the observed mass loss 
in these systems. 

Kinematic modeling of wind-formed, ultraviolet line profiles in CVs suggests mass-loss rates on the order of a few percent 
of the accretion rate (Shlosman, Vitello & Mauche 1997; Knigge & Drew 1997; Prinja & Rosen 1995; Hoare & Drew 1993). 
Guided by this, a value of M Wt total/M a cc(Rdiak) — 0.025 will be adopted for both dynamical pictures below. Note that the 
current generation of numerical and analytical models for line-driven winds is still struggling to achieve such high efficiencies 
(Proga, Drew & Stone 1998; Feldmaier & Shlosman 1999; Feldmaier, Shlosman & Vitello 1999). Centrifugally-driven disk 
winds do not share this pr oble m, but do require the presence of a sufficiently strong, large-scale, ordered magnetic field 
threading the disk (Section [i^J Cannizzo & Pudritz 1988; Pelletier & Pudritz 1992). 

For the mass loss power law index f, values of 0,-1,-2,-5/2 and -3 will be considered for all models. Note that not all of 
these values are equally plausible in the context of a given dynamical picture. This point will be considered in more detail 
in the relevant sections below. Appropriate choices for the remaining outflow parameters (/, r/k and r/ w ) are also made and 
explained in these sections. 

Finally, the following system parameters are adopted throughout: Macc(Rdisk) — 10 _8 MQyr _1 , M* = 0.5Mq and 
R* = 0.0I33i?Q. These are appropriate for a non-magnetic CV in a high state, e.g. a nova-like variable like UX UMa. 

6.1 Radiation-driven disk winds 

As already noted above, radiation-driven disk winds are expected to belong to the I — 1 family of models. The terminal 
velocities in these types of outflows are typically on the order of the (local) escape speeds on the disk (e.g. Proga, Stone & 
Drew 1998; Feldmaier, Shlosman & Vitello 1999), so it is reasonable to take / = V2- 
Two extreme possibilities are considered for the efficiency parameters 7/j, and r^: 

(i) Minimal models [t]b = T)k — 0]: For this choice of parameters, the effect of the outflow on the disk is minimized. The 
minimal model is appropriate if the wind is powered predominantly by an external energy source, such as the boundary 
layer between the disk and the central star. However, somewhat counter-intuitively, it may also be relevant if the disk does 
power the wind. This is because a radiation-driven outflow is quite unlikely to drain energy from the disk's radiation field in 
a manner that would observationally mimic a cooler (blackbody) disk. Most importantly, in a line-driven wind, the transfer 
of momentum and energy from the radiation field to the wind often takes place in a few strong lines. In this case, the overall 
continuous disk spectrum is hardly affected and the minimal model is likely to provide the most relevant benchmark for 
comparisons with observations, fl 

(ii) Maximal models [r)b — T)k — 1] : For this choice of parameters, the effect of the outflow on the disk is maximized. En- 
ergetically, the maximal model corresponds to the assumption that the disk alone is responsible for powering the outflow. 
However, as just noted above, the minimal model may nevertheless be preferable to describe the observational effects of a 
purely disk-driven wind. The maximal model is included mainly to set an upper limit on the effect of a radiation-driven 
wind on the disk's temperature distribution. 

The radial effective temperature distributions corresponding to the minimal and maximal models are shown in Figures hi 
and bJ. In all of these models, the presence of the outflow results in reduced temperatures over the full range of radii, which 
simply reflects the removal of mass and angular momentum: the mass-losing disk needs to exert less viscous torque than a 
conservative to transport material towards the accretor. As a result, the viscous dissipation rates and effective temperatures 
are also lower. 

In most cases, the differences between T e ff(R) with and without mass loss are small everywhere, with 
—AT = — (T e ff — T e ff ;no W i„d) ^ 10 3 K. Such minor deviations are not observable with current instruments and techniques. 

T It is worth noting that, even in a hypothetical "grey" wind (i.e. a radiation-driven outflow with frequency-independent opacity), 
energy is extracted from the radiation field by red-shifting the entire incident spectrum. This is not, in general, equivalent to (or closely 
approximated by) replacing the incident spectrum by one that corresponds to a lower temperature. An excellent discussion of momentum 
and energy transfer in radiation-driven stellar winds is given by Gayley, Owocki &: Cranmer (1995). 
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The only exceptions are the r)t, — T]k = 1, £ < —2 models. For these parameter combinations, the effective temperature 
distribution drops precipitously just beyond the stellar surface. The reason for this behaviour is easy to understand. In these 
models, much or most of the mass lost to the wind is carried away from the inner disk, where the Keplerian and escape 
velocities are highest. Thus the wind binding and kinetic luminosities are largest in these cases, and, for ryt = r\k — 1, must 
be supplied entirely by the disk. As a result, most of the accretion energy that is liberated in the innermost disk regions is 
needed to power the outflow. 

How reasonable are these rjb — T]h — 1, £ < —2 parameter combinations? Intuitively, it seems inevitable that, in a radiation- 
driven accretion disk wind, mass loss should occur preferentially from the hot, luminous, inner disk regions. Dynamical models 
appear to confirm this basic expectation (Proga, Stone & Drew 1998; Feldmaier, Shlosman & Vitello 1999). It is particularly 
worth noting that, in a conservative disk (and, to an equally good approximation, also in the mass-losing disks considered 
here), T e jj oc R~ 3 ^ 4 for R » _R*. This implies that £ = —3 corresponds to the plausible situation in which the mass loss 
rate per unit area is proportional to the rate at which energy is dissipated and radiated locally. In fact, the analytic model 
of Feldmaier, Shlosman & Vitello (1999) predicts rh w (R) oc R~ 29 for R>5Rwd (and rh w (R) oc i? _1 closer to disk centre). 
Thus £ < —2 is a reasonable choice for a line-driven disk wind. 

However, as noted in Section pi the rjb — rjh = 1 maximal model is likely to overestimate the effects of radiation-driven 
mass loss on the disk's continuum spectrum and on the effective temperature distribution that would be inferred for such a 
disk from observations. Moreover, since deviations from the standard case are so highly localized in the innermost disk regions 
in the £ = —3 case, it is highly unlikely that even this model would be observationally distinguishable from a conservative 
disk. Even the most promising technique to infer T(R) from observations, which relies on the inversion of eclipse light curves 
of high-inclination systems (Home 1985), cannot currently provide the high spatial resolution that would be required to make 
such a distinction (A_R « 1 R t ). A final problem is that an optically thick, hot boundary layer between the inner disk and 
the stellar surface could effectively mask this extremely narrow break in the disk's temperature distribution close to the WD 
(Popham & Narayan 1995). 

Thus radiation-driven accretion disk winds are unlikely to leave an observable imprint on the disk's effective temperature 
distribution. 

6.2 Centrifugally-driven disk winds 

The most important parameter that must be fixed to describe a centrifugally-driven outflow within the present formalism 
is the length of the wind's effective lever arm I. As noted in Section kl in this type of outflow IR corresponds to the Alfven 
radius, Ra- 

In the absence of external energy sources, a firm lower bound on I is obtained by noting that I 2 > 3/2 is required 
for such a wind to even overcome its initial binding energy (Section pi). A rough upper bound follows from the fact that a 
centrifugally-driven wind alone - i.e. without any help from viscous torques - can, in principle, drive mass accretion at a rate 
of about (e.g. Cannizzo & Pudritz 1988; Pelletier & Pudritz 1992) 



M acc (RaY 



Mwind \ R J 



(25) 



For M w ,totai I ' Macc(Rdisk) — 0.025, this corresponds to / = 6.3. This is an upper limit, because observations of disk radius 
changes in erupting dwarf nova systems indicate that viscous angular momentum transport must be significant and is probably 
dominant in CVs (e.g. Livio 1999). One arrives at the same conclusion by noting that a purely wind-driven disk, in which no 
angular momentum is transported by viscosity and instead all is extracted by an outflow, would actually be dark. This was 
most recently pointed out by Spruit (1996) and also follows trivially from the fact that there can be no viscous dissipation 
without viscosity: with v — in Equation (hoj) , D(R) and T(R) must also vanish. This upper bound on I is only approximate, 
mainly because boundary effects have been ignored in the derivation of Equation E5l With both of these limits in mind, I = 5 
is adopted in this section to illustrate the effect of a centrifugally-driven wind on the disk. 

Before presenting the results for models of this type, a few preliminary remarks are in order. First, it is of interest to 
consider the strength of the large-scale, poloidal, magnetic field that is required by centrifugally-driven outflows. Following 
Cannizzo & Pudritz (1988), one estimate can be obtained by assuming equipartition between the magnetic and gas pressures 
in the disk photosphere, i.e. 

f^ = pel (26) 

where p is the mass density and c s is the sound speed. Using the standard a-disk solutions (e.g. Frank, King & Raine 1992), 
one finds equipartition field strengths of B eq ~ 2000(-R/10iiVi5) 21/ ' 16 G for the system parameters adopted here. This estimate 
is probably close to an upper limit: even though equipartition of gas and magnetic pressures may be a good approximation 
in the disk atmosphere, it is not necessary for all of the magnetic flux to be in the form of the large-scale field that is needed 
to launch the wind. An estimate of the minimum strength required for this field (denoted B p ) is roughly given by 
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- -ta^F^. (27) 



B p M w i n d H 

~BT q ~ a M acc R 

where H is the disk scale-height and a is the usual Shakura-Sunyaev viscosity parameter (c.f. Pringle 1993). This is larger 
by a factor I compared to Pringle's formula because the wind velocity has been taken to be of order /-times the Keplerian 
velocity. Equation hl\ shows that the large-scale poloidal field can, in principle, be smaller than B eq by factors of 10-100. H 

Second, it should be stressed that the origin of the large-scale magnetic field, as well as the mechanism for maintaining 
it (with a topology suitable for driving an outflow), have not been - and need not be - specified in the present formalism. 
However, it is important to keep in mind that these issues are central to the question of whether centrifugally-driven winds 
can actually exist in real systems (see, for example, Tout & Pringle 1996; Reyes-Ruiz & Stepinski 1996; Lubow, Papaloizou 
& Pringle 1994ab; van Ballegooijen 1989). 

Third, the ratio Ra/R = I might, in general, itself be a function of radius R on the disk (e.g. Pelletier & Pudritz 1992). In 
fact, only £ = —2 is consistent with I = constant in a purely wind-driven, Keplerian disk. This is because angular momentum 
has to be removed from such a disk in just the right way to allow it to accrete while still remaining Keplerian. In "mixed" 
disks, such as those considered here, any combination of / and £ is, in principle, viable, so long as it does not imply the removal 
of more angular momentum than is available. Partly because of this flexibility of mixed disks, the additional complication 
involved in considering, say, a power law form for / did not seem worth pursuing here. If necessary, it would be easy to add 
this, albeit at the cost of at least one additional free parameter. Note that a mixed disk with I — constant and £ = —2 should 
be expected to mimic a conservative disk with a lower accretion rate. 

Fourth and finally, the efficiency parameters r/t and r/k are both set to zero for the models discussed in this section. This 
is because a centrifugally-driven wind already extracts enough energy from the disk to overcome its binding energy and attain 
terminal velocities well in excess of the local escape velocities (i>oo — lv esc ; c.f. Section |3|). 

With these remarks out of the way, the results for this class of models are shown in Figure g. In all cases, the effective 
disk temperatures have been affected more drastically than in their I — 1 counterparts. This was to be expected, since the 
outflows here carry away more angular momentum and hence remove more accretion luminosity in the form of rotational 
kinetic energy. A useful way to compare these two types of models is to note that in both the outflow is accelerated to about 
I times the local escape velocity. Thus a centrifugally-driven wind extracts roughly I 2 times more energy from the disk than 
even a maximal radiation-driven wind (unless the latter also were to attain Woo » v eac ). 

The effects of a centrifugally-driven outflow on T e ff(R) depend on £ in a straightforward manner. Models in which mass 
loss does not decline too rapidly with radius (£ > —2) slightly steepen the slope of the temperature distribution in the outer 
disk region. By contrast, models with £ < —2 tend to flatten the temperature distribution further in and eventually again 
produce a sharp break in the effective temperature distribution near disk center (but note that even for the £ = —3 solution 
1 72. | = 0.027 « 1). For £ = —2 exactly, temperatures are reduced across the disk by an essentially constant factor. 

In dynamical models of centrifugally-driven winds, £ is related to the radial distribution of the large-scale magnetic field 
that threads the disk. It is therefore hard to rule out any particular value for this parameter a priori. The model of Blandford 
& Payne (1982), for example, corresponds to £ = —2. A potential, empirical restriction follows from the fact that the terminal 
velocity of a centrifugally-driven wind is roughly Voo ~ lv esc , where v e3C is the local escape velocity from a streamline footpoint 
on the disk. If a significant fraction of the wind escapes from the inner disk, as in the £ < — 2 models, one might expect to 
see maximum outflow velocities of about /-times the escape velocity from the central object. This is somewhat higher than 
observed in CVs, where the maximum velocities inferred from wind-formed spectral lines are closer to the escape velocity 
itself (e.g. Prinja & Rosen 1995). However, it would be premature to conclude that such models are therefore ruled out in 
CVs. This is because projection effects will tend to reduce the observed maximum velocities. In addition, the outflow might 
be too dilute to produce much line emission or absorption by the time it reaches its terminal speed. Both effects are seen, for 
example, in the kinematic disk wind model developed by Knigge & Drew (1997) for the nova- like variable UX UMa. 

Figure 3 shows that only models with £ < —2 are likely to be observationally distinguishable from a conservative disk. (For 
£ = —2 the net change in the effective temperatures may be large enough, but the temperature distribution closely resembles 
that of a conservative disk with a lower accretion rate.) A CV disk that drives such a centrifugally-driven wind would have 
a redder overall spectrum and would produce non-standard eclipse light curves if observed in a high-inclination system. It 
is interesting to note that such signatures may already have been seen in CVs. Spectroscopically, Long et al. (1994) showed 
that a "missing" or cooler than expected inner disk could remove the discrepancies between Hopkins Ultraviolet Telescope 

8 Note that even if the strength of the large-scale field were close to B eq , its presence would not imply that such a system would be 
classified as a "magnetic CV" . The latter term is reserved for systems in which the disk is disrupted, or prevented from forming, by the 
magnetic field of the WD primary. For reference: a WD dipole field with a surface strength of about 5 X 10 6 G would be required to 
produce a field comparable to B eq at IORwd, (taking the same system parameters as above). A system harbouring such a WD would 
indeed not be a non-magnetic CV. If a disk could form at all, it would be truncated at a radius of about IORwd ■ Even though this is 
not the type of model considered here, the interaction region between a stellar magnetic field and a truncated accretion disks can itself 
be the launch pad for a different type of centrifugally-driven outflow, namely Shu et al.'s (1994) X-celerator wind. 
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observations of the nova-like variable IX Vel and model accretion disk spectra. Moreover, IX Vel is not unique in that model 
disk spectra provide a poor match to its ultraviolet spectrum - standard disk models tend to fail in much the same way 
for other high-state CVs, both nova- likes and dwarf novae in outburst (Wade 1988; Knigge et al. 1997; Knigge et al. 1998). 
Similarly, eclipse maps derived from observations of some high-state CVs, most notably the so-called SW Sex stars, suggest 
that the effective (really: brightness) temperature distribution in the inner disk regions of these systems is considerably flatter 
than expected (Rutten, vanParadijs & Tinbergen 1992). 

Unfortunately, it is far from clear that these observational signatures are really due to the presence of an outflow, and a 
number of alternative possibilities - e.g. magnetic disk truncation, additional radiation sources in the system or unaccounted- 
for radiative transfer effects - are discussed in the articles cited above. Nevertheless, given that CVs (and other disk accretors) 
are known to drive accretion disk winds, attempts to explain discrepancies between standard disk models and observations in 
terms of the effects of these winds certainly have some appeal. 



7 CONCLUSIONS 

The radial distributions of dissipation rate and effective temperature across a Keplerian, steady-state, mass-losing accretion 
disk have been derived and examined. The parametric approach used to achieve this is sufficiently general to be applicable to 
many types of dynamical disk wind models, including both radiation-driven and centrifugally-driven outflows. 

To simplify the calculations, the local wind mass-loss rate was assumed to follow a power law distribution in radius. In 
addition, wind material moving along a streamline originating at radius R in the disk was assumed to effectively co-rotate with 
the disk out to a radial distance IR, so that non-rotating (I — 0), specific angular momentum-conserving (/ = 1) and angular 
momentum-extracting (I > 1) disk winds could all be described within the same framework. To account for cases where the 
wind is partially or entirely powered by energy dissipated within disk, fixed (but arbitrary) fractions of the wind's binding 
and kinetic luminosities were used as cooling terms in deriving the disks effective temperature distribution. In the absence of 
external (non-disk) energy sources, the energetics of the disk+wind system constrain these fractions to be simple functions of 
the length of the lever arm I. With these assumptions and simplifications, the mass-losing disk's effective temperature profile 
could be obtained analytically. 

Effective temperature distributions for CV accretion disks that lose mass to radiation-driven or centrifugally-driven winds 
have been calculated and discussed in detail. For realistic wind mass-loss rates of a few percent centrifugally-driven outflows 
- particularly those in which mass loss is concentrated near disk centre - appear to significantly alter the disk's effective 
temperature distribution. An accretion disk that is affected by the presence of such an outflow could produce spectra and 
eclipse light curves that are noticeably different from those produced by a standard, conservative disk. In fact, existing 
discrepancies between disk models and observations of some systems could probably be reconciled within this framework. 
However, since other mechanisms might also still be able to account for these discrepancies, the latter cannot yet be regarded 
as clear signatures of a disk wind. This uniqueness problem will disappear once a self-consistent and reliable dynamical picture 
of accretion disk winds in CVs and other disk accreting systems has been found. This will allow the free parameters of the 
parametric disk wind model developed here to be fixed. The formalism outlined above can then be used to obtain a simple, 
but reliable estimate of the impact of the wind on the disk's effective temperature distribution. 
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Figure 1. Top panel: Radial effective temperature distributions of disks losing mass to a radiation-driven wind. The outflow is assumed 
to be specific angular momentum-conserving (I = f), with a mass-loss rate of M w total/ ^acc(Rdisk) = 0.025. None of the wind's binding 
or mechanical luminosity is provided by dissipated accretion energy in these "minimal" models (rn = n^ =0). Each line corresponds to 
a different distribution of the local mass-loss rate per unit area with radius, rh m i n d(R) oc R~^, as indicated. Bottom panel: The effective 
temperature differences —AT = — iT^ff — T e j^ ;no m i n d between the mass-losing disk models and a conservative disk. 
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Figure 2. Top panel: Radial effective temperature distributions of disks losing mass to a radiation-driven wind. The outflow is assumed 
to be specific angular momentum-conserving (I = 1), with a mass-loss rate of M w total/ ^acc(Rdisk) = 0.025. All of the wind's binding 
or mechanical luminosity is provided by dissipated accretion energy in these "maximal" models (r/;, = n^ = 1). Each line corresponds to 
a different distribution of the local mass-loss rate per unit area with radius, rh m i n d(R) oc R~^ , as indicated. Bottom panel: The effective 
temperature differences —AT = — T e ^ — T e ^^ ;no m i n d between the mass-losing disk models and a conservative disk. 



© 0000 RAS, MNRAS 000, 000-000 



14 C. Knigge 



Figure 3. Top panel: Radial effective temperature distributions of disks losing mass to a ccntrifugally-drivcn wind. The outflow extracts 
angular momentum from the disk (I = 5) and has a mass-loss rate of M w total/ ^acc{Rdisk) = 0.025. None of the wind's binding or 
mechanical luminosity needs to be provided by dissipated accretion energy in these models (% = n^ =0), since the outflow already 
extracts sufficient energy from the disk to accelerate the ejected material to velocities in excess of the local escape velocity. Each line 
corresponds to a different distribution of the local mass-loss rate per unit area with radius, rh win< j(R) oc R~^, as indicated. Bottom 
panel: The effective temperature differences —AT = — T e ^y — T e ff. no w i na \ between the mass-losing disk models and a conservative 
disk. 
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